Evolution of magnetic field curvature in the Kulsrud- Anderson 

dynamo theory 
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ABSTRACT 

We find that in the kinematic limit the ensemble averaged square of the curvature of magnetic 
field lines is exponentially amplified in time by the turbulent motions in a highly conductive 
plasma. At the same time, the ensemble averaged curvature vector exponentially decays to zero. 
Thus, independently of the initial conditions, the fluctuation field becomes very curved, and the 
curvature vector becomes highly isotropic. 

Subject headings: ISM: magnetic fields — MHD — turbulence — methods: analytical 



It was shown by Kulsrud and Anderson (1992) 
that MHD turbulent dynamo action builds mag- 
netic field energy primarily on scales smaller than 
the smallest turbulent eddy size (which is the vis- 
cosity scale) , but still larger than the resistive scale 
(provided the magnetic Prandtl number is large). 
This result was found under assumption that the 
"kinematic" approximation is valid, i.e. the field is 
weak enough that it does not affect the turbulent 
motions. In this paper we calculate the evolution 
of magnetic field curvature within the framework 
of the Kulsrud- Anderson kinematic dynamo the- 
ory. Such calculations are of great interest because 
a rapid built up of curvature may quickly break 
down the kinematic approximation and change the 
dynamo action on very small scales considerably. 

Following Kulsrud and Anderson (1992) we 
make the following assumptions. We use the 
"kinematic" approximation. We neglect resistiv- 
ity (infinite magnetic Prandtl number limit). We 
assume that the turbulence is incompressible, ho- 
mogeneous, isotropic and static, and we use zero 
correlation time approximation for the turbulent 
motions: 



V a (t,r) = j^ Jv a (t,k)e^d 3 k, (1) 

(2) 



<W,k')V^,k)) = [J k (5 a p - kjtp) 
+ i J/c£ Q7 ^fc 7 ](5(k' — \t)5(t' — t). 



(3) 



(27T) 

{V a (t,k)) = 0, 



Here and below (...) means ensemble average, S a p 
is the Kronecker symbol, £ Q/ 3 7 is the unit anti- 
symmetric tensor, S(t' — t) and <5(k' — k) are the 
Dirac 8- functions, k = k/fc is a unit vector, and we 
always assume summation over repeated indices. 
Functions </& and </& are the normal and the heli- 
cal parts of the turbulence, they depend only on 
the absolute value of k. We make no assumptions 
about the statistics of initial magnetic field. 

Using equations (1) and (3), it is straightfor- 
ward to calculate the correlation tensors between 
velocities and their spatial derivatives, taken at 
the same point of space r but at different time, t 
and t': 

(V a (t',r)V fJ (t,r)} = {r, T l2n)S a pS{t' ~ t), 
(V a Vf3- n ) = ae Q/ 3 7 (5(t / - t), 

{VatpVytf) = (7/5) (BSarySpg ~ <5a/3 7 <5 ) (> (t ' - t), 

(V a V . jS ) = - (V a ,sV , 7 ) , 

(Vc^Vj-st) = ue aiv 5 v p5 T 8(t' - t), (4) 

<y a . Sl V s , Tri ) = (A/15)(7<5 

— <5a/3 7 <5Tij)^(i' — t). 

Here and below, in order to shorten notations, spa- 
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tial derivatives are assumed to be taken with re- 
spect to all indices that are listed after ";" signs, 
e.g. V a -p = dVa/dx/3, V a -,p~t = d 2 V a / dxpOx^ . We 
have also introduced the symmetric tensors 

Sap-yS — SafidyS + S aJ 5ps + SasSpj, 

and the constants (Kulsrud & Anderson 1992) 
4tt 1 



J k 4irk 1 dk, 



3 (2tt) 6 
11 f - 
a = 3(2^)6 J k2J k^k 2 dk, 



1 1 



A 



3 (2tt) 6 
1 1 
15 (2tt) 6 
1 1 



7 (27r) e 



k 2 J k 4:TTk 2 dk, (6) 

fc 4 J k Ank 2 dk, 
k 4 J k Airk 2 dk. 



The equation for the evolution of magnetic field 
B in an incompressible highly conductive fluid 
(Landau & Lifshitz 1983) is 

dB a /dt = V a .fsBf S -Vf 3 B a , !3 . 

We use this formula to derive the equation for the 
evolution of the unit vector b = H/B tangential 
to magnetic field lines: 

db a /dt = Va-pbp - Vfcybabpby - Vpb a -p. (7) 

Now, for a given point of space we assume b is 
known at zero time, t = 0, and we solve this equa- 
tion by iterating it twice in time, similar to the cal- 
culations of Kulsrud and Anderson (1992). Con- 
sidering t > as an expansion parameter, we have 

b(t)=°b + 1 b(t) + 2 b(t), (8) 

where °b is the value of b at zero time, °b = b(0), 
and 1 b(t) oc t and 2 b(i) oc t 2 are the first and 
the second order corrections to b(t), obtained by 
iterating equation (7) twice in time, 

X(t) = I {V a -,p{t")\ ~ WOVV^r 
Jo 

-V ti {t")% a , fi }dt\ (9) 

2 b a (t) = f {V a ., p {t')%{t') - V , 7 (t') x 
Jo 

x [ 1 6 Q (t')°&A + V&/3(0\ 

+ VVM*')] - W) VM*')} df. (io) 



The ensemble averaged curvature 7C = ((bV)b), 
up to second order terms, is 

M*) - (b/3(t)b a , (t)) = (\°b a .,p) 

+ [(\)% a :J i +\( 2 b a ). f3 + (\ 1 b a .J i )] 

= °IC a +t{- (7 7 /5)°/C Q - (4 7 /5)(°6 Q 6 /3;/3 ) 
+ ae Q/ 3 7 /C 7;/3 + (77T/47r) /C Q;/ 3 )3 }. (11) 

Here, we use expansion (8) to find 7C(i). The 
zero order term is °K, a = ( 6/3°6 Q; /3), all first or- 
der terms average out according to formula (2), 
and the second order terms are given in brackets 
[...] on the second line of equation (ll). 1 The fi- 
nal result for JC a {t) was obtained by making use 
of equations (4), (9) and (10). 

Introducing another ensemble averaged vector 
Q = (bdivb), we write the differential equations 
for time evolution of K{t) and Q{t) as 

f = -^/C-^ + a(Vx/C ) + ^A/C, (12) 
at 5 5 47r 

5-- 7 /c-% + «(Vxg) + ^Ac;. (13) 

Ot 5 47T 

Here, the first equation directly follows from for- 
mula (11), and the second equation for the evolu- 
tion of Q can be found by calculations similar to 
the calculations that led to equation (11), i.e. by 
expanding Q{t) up to second order terms, by mak- 
ing use of equations (9), (10) and of equations (4) 
to carry out ensemble averaging. 

In a similar way, and after considerable algebra, 
we find the differential equations describing the 
evolution of the ensemble averaged square of the 
curvature, K 2 = ([(bV)b] 2 ) ^ K 2 , 



dIC 2 

dt 
dD K 



^ 2 + ^Ac + ^AK 2 + ^, (14) 





Air' 



12A 
T 



m -¥^-f D « + S AD - < 15 » 

Here, the differential equations for ensemble av- 
eraged scalars = (div[(bV)b]) = div/C and 



1 Note that we use ( 2 fe Q ;/3) = (^ba).^, which follows from im- 



portant identity ((V c 



a;t3 1 ...fj n V~ f ,S 



t );Tl...Tp) 



0. This 



reflects <5(k' — k) correlation property of the turbulence, 
see eq. (3). Thus, ensemble averaging and taking deriva- 
tives can be exchanged for any second order quantity. 



2 



Dg = (div [bdivb]) = divQ can be found by tak- 
ing the divergence of equations (12) and (13) (see 
the first footnote on page 2). 

We can drop a(V x 7C) and a(V x Q) terms in 
equations (12) and (13) because the helical part of 
the turbulence, Jfe, is usually negligible on scales 
of the smallest turbulent eddy. 2 In this case equa- 
tions (12)-(16) have the following solutions: 





K 


= Qi 


' l " 




Q 


-l 




K 2 ' 




"l" 


D K 


= Qz 













-3ft /S 



3 16 7 t/5 



+ Qa 



+ Q 5 



28 
35 



e -12 7 t/5 + ^ 
47 



4 
5 

8 " 
-19 
19 



,16 7 t/5 



-12 7 t/5 



-3 7 t/5 



Here, the functions Qi(t, r), i — 1,2, 3, 4, 5, are the 
solutions of the same simple diffusion equation 



OQi T)T 



AQi(t,r), 2 = 1,2,3,4,5 (17) 



it=o' 



dt 4tt 
with different initial conditions: 

Qi(0,r) = (l/9)[5/C-4S 
Q 2 (O,r) = (l/9)[/C + 0] t=o . 
Q 3 (0,r) = [/C 2 + (1/133)(48D X - SD g )] 
Q 4 (0,r) = (l/m)[-5D K + ADg\ 
Q 5 (0,r) = (l/63)[D K + D g ] t=0 



t=o> 



Jt=0' 



We see that the ensemble averaged square of 
the curvature, /C 2 , exponentially grows with rate 
I67/3. This is even faster than the rate of mag- 
netic field energy growth, 27 (Kulsrud & Ander- 
son 1992). At the same time, the ensemble av- 
eraged curvature vector, /C, exponentially decays 
with rate —37/5. Therefore, 7C becomes highly 
isotropic. According to diffusion equation (17), 
we find that K? and 7C become homogeneous on 
scales L after a diffusion time ~ A-kL 2 /tjt- 

Let consider the case when the initial mag- 
netic field is constant in space, B = const and 
b = const. Then Qi(t, r) = 0, i = 1, 2, 3, 4, 5, and 



2 E.g. in a galaxy, if the turbulence is Kolmogoroff, J k oc 
k~ 7 and J k oc fc^ 13 / 3 (Kulsrud & Anderson 1992), these 
smallest turbulent scales make the largest contributions to 
a and 7, see eqs. (6). 



we find that K 2 = (3A/4 7 ) (e 16 ^? 5 - l). There- 
fore, even if there is no initial curvature, the field 
quickly becomes very curved. The curvature first 
develops linearly in time, K? w 12 At/ 5, because 
of second order spatial derivatives of the turbu- 
lent velocities, represented by the "battery" term 
12A/5 in equation (14). Second, at time t ~ 5/I67 
the averaged curvature reaches the smallest turbu- 
lent eddy scale size, ^ c ddy ~ 7/ A, and afterwards 
exponentiates rapidly. 

It was recently suggested by Steven Cowley, 
on the basis of numerical simulations, that on 
scales smaller than the smallest turbulent eddy 
the turbulent motions can only stretch magnetic 
field lines. It was argued that the rapid expo- 
nential increase of the fluctuation field wavenum- 
ber k = {k 2 ^ + k 2 ) 1 / 2 found by Kulsrud and An- 
derson (1992) is mainly due to the increase of 
the wave number perpendicular to the mag- 
netic field lines, while the parallel wave number 
fcii stays approximately equal to the smallest eddy 
size, so that k± ^> fc|| . The results we obtained, 
that the ensemble (volume) averaged curvature 
squared, JC 2 , increases exponentially in time, could 
be consistent with the folding nature of small-scale 
fields. For example, the curvature could be com- 
parable to the smallest turbulent eddy size scale 
over the bulk of the volume, and be enormously 
large over the remaining small region of the vol- 
ume (Schekochihin & Cowley 2001). 

The consequence of the rapid development of 
curvature on size scales smaller than the smallest 
turbulent eddy size scale is that the Lorentz force 
and viscosity stress on these very small scales be- 
come important before energy equipartition, be- 
tween the field and the turbulence, occurs on the 
eddy scale. Thus, the magnetic energy at these 
very small scales may be suppressed at least to 
some extent. 
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